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Abstract 

By mapping the Hamiltonians of the two-mode and 2-photon Rabi models to 
differential operators in suitable Hilbert spaces of entire functions, we prove that the 
two models possess entire and normalizable wavefunctions in the Bargmann-Hilbert 
spaces only if the frequency oj and coupling strength g satisfy certain constraints. 

This is in sharp contrast to the quantum Rabi model for which entire wavefunctions 
always exist. For model parameters fulfilling the aforesaid constraints we deter¬ 
mine transcendental equations whose roots give the regular energy eigenvalues of 
the models. Furthermore, we show that for k > 3 the /c-photon Rabi model does 
not possess wavefunctions which are elements of the Bargmann-Hilbert space for all 
non-trivial model parameters. This implies that the k > 3 case is not diagonaliz- 
able, unlike its RWA cousin, the /c-photon Jaynes-Cummings model which can be 
completely diagonalized for all k. 

PACS numbers: 03.65.Ge, 02.30.Ik, 42.50.Pq. 


1 Introduction 

Since Braak’s breakthrough work on the quantum Rabi model [1], there have been renewed 
interest in analytic solutions of this model [2, 3, 4, 5, 6, 7, 8, 9] and its multi-quantum 
and multi-level generalizations [10, 11, 12, 13, 14, 15, 16]. Prior to the works on analytic 
solutions, the fc-photon Rabi model was investigated mainly via numerical diagonalization 
techniques (see e.g. [17, 18]). It was argued in [17] that this model is non-diagonalizable 
and thus ill-defined for k > 2. This is believed to be related to the impossibility of higher 
order squeezing in the traditional sense in that if*a’^ — if{a^)’^, with ^ being complex 
parameter and (a) boson creation (annihilation) operators, is not essentially self-adjoint 
for fc > 2 [19]. For k = 2 (the 2-photon Rabi case), the constraint condition \2g/u:\ < 1, 
with g being the coupling strength and u the frequency, became present [11, 12] if the 
(Bogoliubov) transformations used in solving the model are to make sense. The situation is 
the same in the 2-mode Rabi case, but with the model parameters satisfying the constraint 
\g/oj\ < 1 [11, 12]. These results indicate that the 2-mode and /c-photon Rabi models are 
qualitatively quite different from the quantum Rabi model. 
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The aim of this paper is to provide an independent proof to the following (somewhat 
unexpected) results based on the application of the Bargmann-Hilbert spaces: (i) the 
2-mode and 2-photon Rabi models are dehned in suitable Hilbert spaces of entire analytic 
functions only if uj,g satisfy the constraints \g/uj\ < 1 and \2g/u\ < 1, respectively; 
and (ii) for k > 3 the /c-photon Rabi model is not dehned for all non-trivial model 
parameters to the extent that it does not possess normalizable entire wavefunctions in 
the Bargmann-Hilbert space. We build our proof from the algebraization and partial 
diagonalization of the boson degrees of freedom. For the 2-mode and 2-photon Rabi 
models with model parameters satisfying the aforesaid constraints, we also determine 
transcendental equations whose roots give the regular energy eigenvalues of the models. 


2 Two-mode Rabi model 

The Hamiltonian of the two-mode quantum Rabi model reads 

H = oj{a\ai + 0202 ) + A(T;j -(- g a^ialal + 0102 ), (2.1) 

where (a) are creation (annihilation) operators of boson modes with frequency u (we 
have assumed that the boson modes are degenerate with the same frequency), (JzjCrx are 
Pauli matrices describing two atomic levels separated by energy difference 2A, and g is 
the spin-boson interaction strength. The Hilbert space of the model, T-Lb 0C^, is inhnite 
dimensional, where Tib is the Hilbert space of the boson and is the spin space. 


2.1 Partial diagonalization 


The bosonic fields allow an algebraization. Introduce the operators K±,Ko, 

K-\- = K— = 0i02) A"q = T 02^2 T !)• (2’2) 

These operators form the usual sm( 1, 1) Lie algebra. This algebra has an inhnite-dimensional 
unitary irreducible representation (i.e. the positive discrete series) labeled by the 
Bargmann index k. For the two-mode bosonic realization (2.2) of sm( 1, 1) that we re¬ 
quire here the Bargmann index k can take any positive integers or half-integers, i.e. 

K, = 1/2,1,3/2, • • •. Thus by means of the sm( 1, 1) representation we have decomposed 
the Fock-Hilbert space Tib of the system into the direct sum of inhnite number of sub¬ 
spaces labelled by k = 1/2,1,3/2, • • •. Within the subspace the Hamiltonian is 
given by 

= 2a; (^Ko - + Aa, + gax{K+ + K_). (2.3) 

In other words, the algebraization allows a partial diagonalization of the Hamiltonian 
(2.1) by bringing it into block-diagonal form, 

i/ = 0 R-F), (2.4) 

K 
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where acts in the mutually orthogonal subspaces ( 8 ) with hxed k. Thus the 
problem of diagonalizing the two-mode Rabi model (2.1) is reduced to that of diagonalizing 
each in the corresponding subspace 0 separately. 

The Hamiltonian possess a Z 2 symmetry (parity), PP = where P = 
^nT(Ko-K) 0 jg parity operator in the subspace 0 C^. Thus each "Hj) 0 splits 
into two invariant subspaces 0 |±) labelled by the eigenvalues ±1 of P. This parity 
invariance can be used to partially diagonalize This is seen as follows. Dehne the 

unitary operator 


— i ^ ^ 

V2\ T -T 


rp _ ^ilz{Ko-K,) 


(2.5) 


Working in a representation dehned by ax diagonal, we have 


t/t fjM u 


0 \ 

0 j ’ 


( 2 . 6 ) 


where for fixed k 

//i") = 2a; {k, - 1) + g{K+ + K.) ± AT. (2.7) 

act in two mutually orthogonal subspaces 0 |±) with hxed parity. 


2.2 Wavefunctions and constraints for uj and g 

In the same way as the differential realization of boson operators in a Hilbert space of 
entire analytic functions of growth (|, 1 ) [ 20 ], we can represent the continuous boson 
degree of freedom {K± o) as differential operators in a Hilbert space of entire functions of 
growth ( 1 , 1 ) [ 21 ]. 

Let Bk be the Hilbert space associated with k whose elements are entire analytic 
functions, called Bargmann-Hilbert space Its inner product is dehned by [21] 

{f, 9 )K = f f{z)g{z)dg^{z), dg^{z) =-\z\^^-^Ki^^{2\z\) dxdy, ( 2 . 8 ) 
J TT 2 

where dg^^i^z) is the measure and K^{z) is the modihed Bessel function of the third kind 
which has the Mellin transform, = r(s-|-2Q;)r(s-|-2/3). Then 

f{z) belongs to B^ if and only if ||/||^ = (/, /)« < 00 . It is now not difficult to see that 
if f{z) = En=o CnZ"^, then 

OD 

ll/ll^ = E Iwpn! (n + 2K-1)!. (2.9) 

n=0 

Every set of coefficients Cn for which the sum on the right hand side converges dehnes 
an entire function f{z) G An orthonormal set of basis vectors in B^ is given by 
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the monomials /\Jn\{n + 2 k — 1)!|. In this basis K±,Ko (2.2) are realized as the 

single-variable differential operators [11, 21] 

Kq = z-^ + K, K+ = z, K_ = z-^+ 2 k-^. (2.10) 

(JjZ LLZ (ajZ 

The operator T can be realized as T = which acts on elements f{z) of as 

{Tf){z)=f{-z). 

By means of the differential representation (2.10), we can express the Hamiltonian 
(2.7) as the differential operator in 

The corresponding time-independent Schrodinger equations are 

+ 2{uz -1- ± + gz + 2u (^K - -E±| = 0. (2.12) 

Here we have written E± since in general the spectra of are not the same. As in 
the Rabi case [22], E± belong to the spectra of if and only if for these values of E± 
solutions to the above differential equations are entire functions belonging to In other 
words, we are seeking solutions of the form 

CXD 

= E z". (2.13) 

n=0 

which converge in the entire complex plane and are elements of 

Substituting (2.13) into (2.12), we obtain the 3-step recurrence relation, 

sf+ c„±sf = 0, 

sy, + cjsj + cjsy, = 0, n>l, (2.14) 



(2.15) 

(2.16) 


when n ^ oo. Thus the asymptotic structure of solutions to the n > 1 part of (2.14) de¬ 
pends on the Newton-Puiseux diagram formed with the points Po(0, 0), Pi(l, —1), ^ 2 ( 2 , —2) 
[23]. The characteristic equation of the 3-term recurrence relation is given by t^-|-^t-|-l = 
0, which has two solutions ~ ^ \/^ ~ Thus we have two cases to consider. 
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Case (1) ; 1^1 < 1. In this case, we have two distinct real roots ti = ^ ["—1 + \Jl — {g/ujY 
t 2 = —^ 1 + — {g/uY , and |ti| < \t 2 \- The Perron-Kreuser theorem (i.e. Theorem 


2.3 of [23 


gives the asymptotic behaviour of two linearly independent solutions 


1. ■-^n+l,r , —1 

lim ^ trU , 


n—)-oo 


r = 1,2. 


(2.17) 


So solution is minimal while S ^2 is dominant. From (2.9), for an entire solution in 
the sum 

CXD 

^ \S^\^nlin + 2K-l)\ (2.18) 


n=0 


must converge. Using the ratio test, 

|^^+iP(n + l)!(n + 2/t)! 
\S^\'^ n\ {n + 2k — 1)\ 


= \tr 


(2.19) 


It is easily seen that |t 2 p > 1- We can show that |tip < 1. Indeed, if we had |tip > 1, then 
we would end up with — \g/u)\ > ^Jl + \g/u\ which is impossible for the non-trivial 
case S' 7 ^ 0. 

It follows that the sum (2.18) converges for the minimal solution and thus the 
corresponding wavefunctions (2.13) are elements of 13^. 

Case (2); 


and |ti| = |t 2 


> 1. In this case, the two roots ti 2 are complex conjugate to each other 


= 1. Applying the Perron-Kreuser theorem, we have 


lim sup (I S' 

n^oo V ' 


n\ 


= 1 


( 2 . 20 ) 


for all non-trivial solutions of the 2nd equation of (2.14). Thus given e > 0, there exists 
N{e) G N and an inhnite set I of indices ^ > N{e) such that > 1 — e, i.e. 


> (1 - eY/^\. So we have 


Y, \SY\^ n\ {n + 2 k - 1)\ > Y \SY\^i\ {i + 2K-l)\ 

n=0 e&I 


> H (1 - e’ 

eel 


21 {^ + 2K-l)\ 

' J\ ' 


Noting that when £ —)■ cx), e —)■ 0, we have 


1 - 71 {^ + 2K — 1 )! 


( 2 . 21 ) 


( 2 . 22 ) 


This means the series on the right hand side of (2.21) diverges. By the comparison test, 
the series on the left hand side of (2.21), i.e. the sum (2.18) diverges for all non-trivial 


solutions of the 3-term recurrence relations. Thus when 
has no entire wavefunctions which belong to 13^. 


> 1 the two-mode Rabi model 
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2.3 Energy spectrum 


We now proceed to find regular energy eigenvalues E± corresponding to \g/oj\ < 1 and 
the minimal solutions (and thus to the entire wavefunctions (t)±{z)). We 

follow a procedure presented in [24, 2, 12, 14] that uses the relationship between minimal 
solutions and inhnite continued fractions [23]. 

By the Pincherle theorem (i.e. Theorem 1.1 of [23]), the ratios of successive elements 
of the minimal solutions can be expressed as continued fractions. 


C=bmm 
Q± ^ ^n+1 


-^n+l -^n+2 

C di 

n+l“ ^n+2“ ^n+3“ 


(2.23) 


which for n = 0 reduces to 


*50^ = 




ibmm 




Cibmm 

*^0 


/^ib 

Oi — O2 — <^3 “ 


■y± 


Dl 


'± 


(2.24) 


The ratios = 


a__ 

-izhmin 


involve although the above continued fraction expressions 


are obtained from the 2nd equation of (2.14), i.e the recurrence (2.14) for n > 1. On 
the other hand, for single-ended sequences such as those appearing in the inhnite series 


expansions (2.13), the ratios Sq = of the hrst two terms of a minimal solution are 

unambiguously hxed by the n = 0 part of the recurrence (2.14), that is. 


= -Ct 


1 

2gn 


i?± =F A - 2a; 



(2.25) 


In general, (2.24) and (2.25) can not be both satished i.e. the 5^ computed from (2.24) 
are not the same as that from (2.25) for arbitrary values of the recurrence coefficients Cn 
and Dn- Thus general solutions to the recurrence (2.14) are dominant and are usually 
generated by simple forward recursion from a given value of Sq . Physical meaningful 
solutions are those that are entire and normalizable with respect to the the Bargmann- 
Hilbert space norm given in (2.8). They can be obtained if E± can be adjusted so that 
equations (2.24) and (2.25) are both satished. Then the resulting solution sequences 
S^{E±) will be purely minimal and the corresponding power series expansions (2.13) will 
converge in the whole complex plane and be elements of Equating the right hand sides 
of (2.24) and (2.25) yields implicit continued fraction equations for the regular spectrum 


E. 


±, 


0 - Co 


Di 


Dl 


(2.26) 


^2 “ ^3 “ 

Here CJ.BJ are dehned as functions of E^ in (2.15), and thus (2.26) contain E± as 
parameters. These are transcendental equations for the determination of the regular 
energies E± of the two-mode Rabi model. Only for the denumerable inhnite values of E± 
which are the roots of (2.26), do the power series (2.13) give convergent and normalizable 
solutions to the diherential equations (2.12). 
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3 /c-photon Rabi model 


The fc-photon Rabi model is a natural generalization of the quantum Rabi model. Its 
Hamiltonian is given by 


H 


(jja^a 


A (T, 


9 


)^ + a^ 


(3.1) 


where k = 1, 2, • • • is a positive integer, (a) are creation (annihilation) operators of 
a boson mode with frequency uj, (7 z,(Tx are Pauli matrices describing two atomic levels 
separated by energy difference 2A, and g is the coupling strength. As in the 2-mode case, 
we will let denote the Hilbert space of the /c-photon Rabi model, where Tib is the 

Hilbert space of boson number states and is the spin space. The k = 1 case of (3.1) 
gives the Hamiltonian of the Rabi model. 


3.1 Partial diagonalization 

The boson helds allow an algebraization. Introduce three operators 



These operators form a polynomial algebra given in [25]. The polynomial algebra has an 
inhnite-dimensional unitary irreducible representation labelled by the parameter q, which 
for k = 2 reduces to the well-known positive discrete series of sm( 1, 1). For the single-mode 
bosonic realization (3.2), the parameter q takes the k allowed values [25], 

1 k -\- \ 2k -|- 1 [k — l)k -|- 1 ^“5 

Thus by the algebraization of the bosonic held, we can decompose the Hilbert space Tib 
into a direct sum of k independent subspaces T-Ll labelled by q. Within the subspace 
with hxed g, where n = 0,1, 2, • • •, the Hamiltonian is given by 

H^<i) = ku + A -1- (Q+ -h Q-) . (3.4) 

In other words, the algebraization allows a partial diagonalization of the Hamiltonian 
(3.1) by bringing it into block-diagonal form, 

// = 0 //<», (3.5) 

where Ffacts in k mutually orthogonal subspaces with hxed q. Thus the problem 

of diagonalizing the /c-photon Rabi model (3.1) is reduced to that of diagonalizing each 
//(<?) 

in the corresponding subspace T-Ll ® separately. 

The Hamiltonian possess a Z 2 symmetry (parity), Pq = where Pq = 

<8) (Jz is the parity operator in the subspace Pi ® Thus each Pi 0 splits 
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into two invariant subspaces T-Ll ® |±) labelled by the eigenvalues ±1 of Pq. This parity 
invariance can be used to partially diagonalize This is seen as follows. Dehne the 

unitary operator 


U = 


1 

71 


1 1 
T -T 


T = 


gi7r(Qo-g) 


(3,6) 


Working in a representation dehned by diagonal, we have 


U 


i/7 0 \ 

0 iii®^ ) ’ 


(3.7) 


where for fixed q 

ii7 = ku (Qo -^)+ 9^ {Q+ + Q-)±AT (3.8) 

act in two mutually orthogonal subspaces T-Ll ® |±) with hxed parity. 


3.2 Bargmann-Hilbert space and 3-term recurrence relations 

The continuous boson degree of freedom Qi.o of the /c-photon Rabi model can be realized 
as differential operators in a Hilbert space of entire functions of growth (1,1) [26]. Let 
Bq denote the Hilbert space associated with index g, called Bargmann-Hilbert space Bq. 
Similar to the two-mode case, the inner product in Bq is dehned by 

{f,9)q = f f{z)g{z)dnq{z), dfiq{z) = dxdy. (3.9) 

J 71 

f belongs to Bq if and only if ||/||g = {f,f)q < oo. Now let f{z) be an entire function 
with the power series f{z) = Then it can be easily checked using r(s) = 

/o°° for Re(s) > 0 that in terms of the expansion coefficients 

OO 

ll/llg = E \cn\^[k{n + q-l/k^)]\ (3.10) 

n=0 

Every set of coefficients c„ for which the sum on the right hand side converges dehnes an 
entire analytic function f E Bq. An orthonormal set of basis vectors in Bq is given by the 
monomials \^z^ j 7[/c(n + q — l//i;^)]!|. In this basis, Q±fi are represented as single-variable 
differential operators [26], 



Moreover the operator T can be realized as T = which acts on elements f{z) of 

Bq as {Tf){z) = f{-z). 



Using the differential realization we can equivalently write (3.8) as the single-variable 
differential operator in Bg, 


= kuiz-^ + q-^]±A + gz 

\ dz J 

/ d {j-l)k + V 

+gk'^z ^Y[[z— + q-- 


i=i 


dz 


k^ 


(3.12) 


The corresponding time-independent Schrodinger equations are 


gk’^z-^li 

i=i 


^ d {j-l)k + l\ 

-P- 


+kzj (^zT + 9 - pj ± Ae“"E - B±| i>^(z) = 0. 

r(i) 


(3.13) 


Here we have written E± since in general the spectra of 11^ are not the same. Solutions to 
the differential equations must be analytic in the whole complex plane and normalizable 
with respect to the norm given in (3.9) if belong to the spectra of In other 

words, we seek solutions of the form 


n=0 

which converge in the entire complex plane and are elements of Bg. 

Substituting (3.14) into (3.13), we obtain the 3-step recurrence relation, 

At + .yA'„±=0, 

A'y, = 0, n> 1, 


(3.14) 


(3.15) 


where 


/i± — 


HT = 


±(-1)’" A - E± +ku (^n +q - 
gk'^ n-=i (n + 1 + g - 


k^ n-=i (n+l + g-^^4^ 
The coefficients A^, have the behavior 


~ an", r-^bn^ 


when n ^ oo, where 


u 


a = 


gk^ 


- 1 ’ 


a = —k -|- 1, b = —fd = —k. 
fz 


(3.16) 


(3.17) 


(3.18) 


Thus the asymptotic structure of solutions to the n > 1 part of (3.15) depends on the 
Newton-Puiseux diagram formed with the points Po(0, 0), Pi(l, —k + 1), ^ 2 ( 2 , —k) [23]. 

For /c = 1, (3.15)-(3.18) reduce to the corresponding relations obtained in [5] for 
the Rabi model. It was shown that in this case entire wavefunctions in an appropriate 
Bargmann-Hilbert space exist for all model parameters. 
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3.3 Wavefunctions and constraints for uj and g in the k = 2 case 

For k = 2 the characteristic equation of the 3-term recurrence relation is given by -|- 
-|- I = 0, which has two solutions ti ^2 = ^ Similar to the two-model 

Rabi case, we now have two cases to consider. 


Case (i); 


2 £ 


< 1. In this case, we have two distinct real roots ti = ^ 


-1 -F - {2g/uy 


i — ^ 

^2 - -Tg 

2.3 of [23 


1 -f — {2g/u)^ , and |ti| < \t 2 \- The Perron-Kreuser theorem (i.e. Theorem 
) gives the asymptotic behaviour of two linearly independent solutions 


lim 

n^oo 

n,r 


tr n 


-1 


r = 1,2. 


(3.19) 


So is the minimal solution and K ^2 is dominant. From (3.10), for an entire solution 


in Bq the sum 


y |Atiq2(n + «-l/4)]! 

n=0 


must converge. Using the ratio test. 


A-J+,l"|2(n+l + 9-l/4)|! 


lim 


— 4 \t^ 


(3.20) 


(3.21) 


\K^\‘^[2{n + q-l/4)]\ 

It is easily seen that 4 |t 2 p > 1- We can show that 4|ti|2 < 1 . Indeed, if we had assumed 
4 1 tip > 1 , then we would end up with — \2g/u\ > \J^ + \2g/uj\ which is impossible 
for the non-trivial case 7 ^ 0 . 

It follows that the sum l-^n P [2(n-|-g —1/4)]! converges for the minimal solution 

Thus the corresponding wavefunctions given by (3.14) are entire functions in Bq. 

Case (ii); ^ > 1 . In this case, the two roots ti ^2 are complex conjugate to each other 
and pi I = P 2 I = Applying the Perron-Kreuser theorem, we have 


lim sup (I A'jfl nil" = - 

n^oo i ft I y 2 


(3.22) 


for all non-trivial solutions of the 2nd equation of (3.15). It follows that for a given 
e > 0, there exists N{e) G N and an infinite set of I of indices i > N{e) such that 


(|AtK!) 


1 

" > 4 

2 


e, i.e. |iFp| > Thus we have 


E |App[2(n + g-l/4)]! > ^ |iFpp [2(£ + g - 1/4)]! 
n=o eel 


eel 


(iiy 


Noting that when £ ^ cx), e —)■ 0, we have 


lim - - e 

€^00 \2 


2i 


[2(£ + l + g-l/4)]! 

(U )2 


7^ 0 , 


(3.23) 


(3.24) 
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which means the series on the right hand side of (3.23) diverges. Thus by comparison test, 
the sum J2^=o l-^n P [2(n + g — 1/4)]! diverges for all non-trivial solutions of the 3-term 
recurrence relations and the 2-photon Rabi model has no entire wavefunctions belonging 


to Bq when 




> 1 . 


3.4 Energy spectrum in the k = 2 case 


As in our previous discussion for the regular energy spectrum of the 2-mode Rabi model, 
the coefficients will be minimal solutions = K^i iff they satisfy the continued 

fraction equations 


T^±min pit p± pit 

J^n+1 _ -Pn+1 -Pn+2 ^n+3 

l/'±min 4± _ 4± _ 4± _ 

^n+1 ^n+2 ^n+3 


(3.25) 


which in turn will require that E± be the roots of 


0 = A± 


Bf Bt Bj 
At- At- 


(3.26) 


Here At, Bt are given as functions of E± in (3.26) with k = 2. These are transcendental 
equations whose solutions determine the regular energies E± of the 2-photon Rabi model. 
Only for the denumerable infinite values of E± which are the roots of (3.26), do we get 
entire wavefunction solutions in Bq to the k = 2 version of the differential equations (3.13). 


3.5 /c-photon Rabi model for the k > 3 case 

We now focus on the k > 3 case. Let a be the slope of PqPi and r the slope of P 1 P 2 so 
that a = a and t = /3 — a. Then a = —k -|- 1, r = —1, and thus point Pi lies below the 
line segment PqPi in the Newton-Puiseux diagram. Applying the Perron-Kreuser theorem 
(i.e. Theorem 2.3 of [23]), we have 

mn sup (l/qi («!)?)"= ^ (3,27) 

for all non-trivial solutions of the 2nd equation of (3.15). To have entire solutions which 
are elements of Bq, the sum on the right hand side of (3.10) must converge. To check if 
this is the case, we first note that similar to the k = 2 case, for a given e > 0, there exists 
N{e) G N and an infinite set of / of indices ^ > N{e) such that (\Pe^\ (^!)^) ^ > — e, 

i.e. \Kt\ > So we have 

OD 

ft \KtmHn + q-l/k^)V- > Y. 

n=o eel 

> E 

eel 


\Kft[k{^ + q-l/k^)]\ 


2e 




— e 


[kii + q-l/k"^)]'. 
(H)^ ■ 


(3.28) 
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Now 


lim 

i^OO 


1 


(3.29) 


21 




— e 


[k{i+l + q-l/k^)]\ 

(£!)^ ^ 


This means the series on the right hand side of (3.28) diverges. It follows that the sum 
on the left hand side of (3.28) diverges for all non-trivial solutions of (3.15) and the 
Schrodinger equations (3.13) with k > 3 have no solutions belonging to Bq. 

We may thus conclude that for /c > 3 the /c-photon Rabi model does not have eigen¬ 
functions which are elements of Bq because they are not normalizable with respect to 
the norm given in (3.9). In other words, for k > 3 the Hamiltonian (3.1) can not be 
completely diagonalized in the Hilbert space "Hf, <8) due to the non-normalizability of 
its eigenstates. This is in sharp contrast to the /c-photon Jaynes-Cummings model which 
can be completely diagonalized for all k [27]. The impossibility conclusion for the k > 3 
case agrees with that reached in [17, 18, 19] by using different approaches. 


4 Conclusions 

We have examined the 2-mode and /c-photon Rabi models based on the application of 
algebraizations and Bargmann-Hilbert spaces. We have seen that the algebraization and 
parity invariance allow us to decompose the total Hilbert spaces into direct sums of 
independent subspaces, thus partially diagonalizing the Hamiltonians of the models by 
bringing them into block-diagonal forms. The block-diagonal sectors can be realized as dif¬ 
ferential operators in the Bargmann-Hilbert spaces. We have investigated the eigenvalues 
and eigenfunctions of the block-diagonal sectors by applying the theory of Bargmann- 
Hilbert spaces. We have derived constraints for the frequency u and coupling g for the 
2-mode and 2-photon Rabi models to be dehned in the Hilbert spaces of entire analytic 
functions. We have determined the corresponding transcendental equations whose roots 
give the regular energy spectra of the models. Furthermore we have shown that the k- 
photon Rabi model with k > 3 does not have normalizable eigenfunctions with respect to 
the Bargmann-Hilbert space norm and thus can not be completely diagonalized. 

As shown in the Appendix, the wavefunction expansion coefficients for the 2-mode 
and 2-photon Rabi models are related to orthogonal polynomials. Thus it is expected 
that the regular energies of the two models can be determined as the polynomial zeros by 
a procedure similar to that in [5, 6]. Work in this direction is in progress and results will 
be reported elsewhere. 
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A Orthogonal polynomials 


In this appendix, we show that the expansion coefficients in 4>^{z) (2.13) and (3.14) 

are related to orthogonal polynomials of inhnite degree. 

Let us recall the well-known theorem [28] on relationship between 3-term recurrence 
relations and orthogonal polynomials. It states: the necessary and sufficient condition 
for a family of polynomials {P„(a;)} (with degree = n) in parameter x to form an 
orthogonal polynomial systems is that these polynomials satisfy the 3-term recurrence 
relation 

Pn{x) = {PnX - an)Pn-l{x) - \nPn-2{x), U > I, (A.l) 


where the coefficients I3n, an and are independent of x, f3n Q and A„ 7 ^ 0 for n > 
1. Then {Pn{x)} forms an orthogonal set of polynomials with respect to some weight 
function. 

For the two-mode Rabi case, the expansion coefficients S^{E±) dehned by the 3-term 
recurrence relations (2.14) are related to orthogonal polynomials of inhnite degree in 
energy parameters E±. To see this, dehne P^{E±) by 




P^(E±) 
n\ {n + 2K — 1)! 


(A. 2 ) 


Then in terms of Pj^, (2.14) becomes 


p± _ 

n+l ~ 


E± =F (-1)"'A - 2uj (n + K- 


P^ -n(n + 2K- l)P^_i. 


(A.3) 


It then follows from the above mentioned theorem that P^{E±) are orthogonal polyno¬ 
mials in E± with degree n. Thus (t)^{z) (2.13) are the generating functions for P)^(P), 




pHe^) 


z = 


00 

{ri + 2K-l)\ 


(A.4) 


However, (l)^{z) are entire only if P^{E) dehned in (A.2) correspond to the minimal 
solutions Sn^i{E). 

Similarly for the /c-photon Rabi case, the expansion coefficients K^{E±) dehned by the 
3-step recurrence relation (3.15) are related to orthogonal polynomials of inhnite degree 
in energy parameters E±. This is seen as follows. Dehne P^{E±) by 


K^E^) 



P^{E±) _ 

k-T(n+l + q-^i^^) 


(A.5) 
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where r(a;) is the gamma function in x. In terms of P^{E±), (3.15) become 


P 


n+1 


E± T (-1)”A -kuin + q- — 






-n 

j=i 


k in + q — 


{j - l)fc + 1 ^ 
k^ 


p± 


(A.6) 


Thus P^{E±) defined by (A.5) form a set of orthogonal polynomials and the solutions 
(3.14) are generating functions for P^(P±), 




OO 

E 

n= 


P^iE) 


To nLi k^T [n + l + q 


fc2 ) 


(A.7) 


However ip^iz) given above have hnite radius of convergence (i.e. they are not entire) for 
k > 3. So we focus on the k = 2 case in which entire wavefunctions exist. Note that in 
this case the denominator of (A.5) equals 2"'r(n + l + g — 1/4) x2"'r(n + l + g — 3/4) = 
[2(n + q — 1/4)]!. So the solutions for /c = 2 can be expressed as 




^ Pt{E±) 

i) [2(n + g-l/4)]! 


(A.8) 


When P^{E±) correspond to the minimal solutions K^i{E±), ijj'^{z) are entire functions 
and thus give the wave functions of the 2-photon Rabi model. 
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